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Abstract 

We present here results of numerical calculations demonstrating that the very 
long periodic variability detected in the Cataclysmic Variable FS Aur can be the 
result of eccentricity modulation of a close binary (CB) orbit induced by the pres- 
ence of a third body on a circumbinary orbit. A third component with a substellar 
mass on a circular, relatively close orbit, modulates the mass transfer rate of the 
binary on much longer time scales than periods within the triple system. 
We also report here preliminary results of X-ray observations of FS Aur, provid- 
ing further evidence that it contains magnetic and freely precessing white dwarf. 
These two findings allow us to incorporate the new and previously stressed hy- 
pothesis on the nature of FS Aur into one consistent model. 

1. Introduction 

FS Aur is a short period Dwarf Nova (DN) renown for the presence of a variety 
of uncommon and largely incomprehensible periodic variabilities of brightness and 
radial velocity. The orbital period (OP) is determined from the radial velocity 
(RV) variation of the central peak of emission lines, also referred to as S-wave. 
The widely used method of making RV measurements using the doublc-gaussian 
method reveals a different picture in the case of FS Aur because the wings of the 
lines are dominated by a periodicity that is longer than the OP, which we call the 
long spectroscopic period (LSP). The light curve (LC) is even more complicated. 
The most persistent and dominant LC modulation is the so-called long photomet- 
ric period (LPP). There is, however, a definite relation between OP, LSP & LPP: 
1/Plsp = 1/-Ppp — l/^ippi so tne LSP and LPP seem to emerge from the same 
source within the system. The LC is also highly variable in its appearance, only 
sometimes showing the OP. There is also some non-periodic, stochastic variability 
contaminating the LC on time scales of 1 — 2 x 10 3 sec, much longer than the flick- 
ering, which is also prominent in the LC. The amplitudes of the variabilities range 
from < 0.1 mag flickering to « 0.5 mag for LPP and stochastic variability. How- 
ever, the long term monitoring of FS Aur reveals an additional, discordant, and 
very long photometric period Pvlpp ~ 900 days with an even greater amplitude. 
In FigfTjup to five 2 m dips can be detected in the LC of FS Aur publicly available 
from AAVSO () since it became closely monitored during the last decade. 

The optical mix of a DN is dominated by the accretion disk, and the stellar 
contribution is small or negligible. Thus, the variability of a DN is mostly the 
result of internal structure or conditions in the disk or as a result of a variable 
mass transfer rate. We demonstrate here that the variability with VLPP can 
be induced by a long-term eccentricity modulation of CB orbit, influencing the 
mass transfer rate by a third, sub-stellar component orbiting the binary on a close 
circular orbit with a period much shorter than the observed VLPP. On the other 
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Fig. 1. Long-term light curve of FS 
Aur from the AAVSO data base. 
The open (cyan) points correspond 
to outburst that were duly sepa- 
rated using a V> 15 mag thresh- 
old. That left some of the outbursts 
occurring from deep minima unfil- 
tered, but these events are rare and 
do not affect the period analysis at 
very low frequencies. The time in- 
dicated in the lower axis is in JD 
and in the upper axis is in years. 




;. 2. The upper panel shows the VLPP 
modulation in binary eccentricity as a 
function of the perturber mass. The 
curves correspond to different P?,IPn- 
The thick horizontal line shows the 
VLPP value. Solutions that cross this 
line can explain the VLPP. The middle 
panel shows the perturber mass and or- 
bital radius combinations that result in 
a long-term modulation of the binary 
orbit equal to the VLPP. The lower 
panel shows the amplitude of the bi- 
nary eccentricity perturbation. 



hand, we will show that the system probably contains a freely processing magnetic 
white dwarf (WD) and that LPP and LSP then can form in a coupling region on 
the inner edge of the accretion diskpQ. 

2. FS Aur as a triple system 

The LC of FS Aur is very complex with large nightly variations sometimes 
exceeding 0.5 mag. Even with such a large spread, there are also obvious 2 mag- 
nitude dips in the average brightness of the system. The dips are sharp and the 
variability is not sinusoidal, but the succession of minima appears to be periodic 
(Fig[L]). The period determined by the Discrete Fourier Transform (DFT) pro- 
cedure returns a value around 875 days. With only five cycles observed during 
the last 10 years, during which time we had systematic and relatively dense cov- 
erage, the period determination is very rough. We applied very primitive filters 
for the observational data, not caring about the precision of the period determi- 
nation. The irregular outbursts of FS Aur were excluded by deleting points with 
magnitudes brighter than 15th. Only CCD based observations were taken into 
account and these were predominantly taken in V-band. The unfiltcred data are 
included too, since the color difference is negligible compared to the magnitude 
of the variability The long nightly observational runs were averaged in the final 
LC presentation. The period analysis results were not sensitive to whether the 
nightly data were averaged or not. In hierarchical triple systems, a third light- 
weight body can produce perturbations on a central binary, whose period is long 
compared with the orbital period of the perturber ([2], [3], [4]). A third companion 
prevents the complete circularization of the orbit of a binary system by inducing 
a long-term eccentricity modulation. This long-term modulation is produced by 
the time- varying tidal force of the perturber upon the binary. We performed nu- 
merical analysis to explore the range of possible solutions. The integrations were 
performed using the high-order Runge-Kutta-Nystrom RKN 12(10) 17M integra- 
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tor of [5] for the full equations of motion in the barycenter frame. At each time 
step, the instantaneous eccentricity of the binary is computed using the oscillating 
ellipse approximation (Eq. 2.135 of 0). The basic setup is a binary formed by 
two point masses initially in circular orbit. We used the masses of the binary 
components common for a 86 min CV (Mi = O.7M & M 2 = 0.14M Q ), with the 
orbital period Pi equal to the OP. A third point mass (perturber) moves initially 
on its own circular orbit, farther away and in the same plane. Its mass M3 and 
orbital period P3 are changed across an ensemble of numerical experiments. 

The top panel of Figure[2~1 shows the resulting periods of the long-term modu- 
lation of the binary eccentricity as a function of the mass of the perturber. The 
solid horizontal line corresponds to the VLPP value. For a perturber whose orbital 
period is less than 13.4 binary periods, no solution is possible, as their respective 
curves do not reach up to the VLPP value. For a perturber with periods longer 
than that but shorter than 23 orbital periods, two solutions are possible: one at 
low mass and another at an increasingly larger mass. Finally, a perturber with 
longer periods produces only one solution at the large mass end of the range. The 
curve in the middle panel of Figure[2Dpresents the perturber orbit semi-major axis 
for the solutions that result in the binary system eccentricity modulation of pe- 
riod equal to the observed VLPP value, i.e., solutions which cross the solid line on 
the upper panel. The lower panel of Figure[2~1 shows the maximum amplitude of 
the eccentricity perturbation that can be reached for the corresponding solutions 
presented in the middle panel of this figure. The maximum amplitude is achieved 
for systems which harbor a third body with M3 = 48Mj. 

The solutions depend weakly on the exact parameters of the CB, and a wide 
range of third body masses and orbits will induce the long period effect on the 
CB eccentricity. More precisely, a third body with P3 / Penary period ratios from 
12 to 48, and mass range from a few to a hundred Jupiter masses would be able 
to produce eccentricity modulation of order of 900 days. However, the binary 
eccentricity would be modulated on three different time scales: on the binary 
period Penary ' on the period of third body P3 and on the VLPP. The amplitude 
of each of those modulations would depend primarily on the perturber's mass. 
The LC folded with 875 day period is presented in the top panel of Figl3~Tlcft'). 
In the bottom panel of Fig|3~Tright). the modulation of the binary eccentricity 
is presented for three distinct cases: low, intermediate and high perturber mass 
(left to right panels). The amplitude of the low frequency modulation Ae;, which 
basically corresponds to the VLPP, is measured by the maximum to minimum 
difference of the outline of the corresponding curves. The amplitude of the high 
frequency modulation Ae^, is determined from the width of the strip forming 
the curves. The measure of relative amplitude is expressed as a ratio of these 
amplitudes (© = Ae;/Ae^). The upper panel of Figl3~Tright) demonstrates the 
continuous dependence of the amplitude ratio, G, on the mass of the perturber. 
The three particular cases shown at the bottom of the figure are marked on the 
curve. It is apparent that the amplitude of the low frequency eccentricity variation 
is the strongest compared to that of the high frequency, when the perturber's mass 
hovers around 50 Jupiter masses. Formally, the amplitude of the low frequency 
modulation becomes dominant (0 > 1.0) in the range of masses for the third 
body of 25Mj < .A/3 < 65A/j. The higher mass solutions are also excluded by the 
absence of any IR excess in the spectral energy distribution of FS Aur. The long 
term eccentricity modulation curve produced by a third body in such mass range is 
remarkably similar to the observed LC of FS Aur as can be seen on the lower panel 
of Figj3jleft). Of course, the physics of translating the eccentricity variation to the 
brightness of the accretion disc is not that simple, but the eccentricity modulation 
directly changes the Li distance from the center of the mass of the secondary and 
therefore the mass transfer rate. The eccentricity modulations on much shorter 
time scales have smaller amplitudes, but may be reason for the erratic behavior 
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Fig. 3. left) In the upper panel, the relative amplitude of the low frequency to the 
high frequency ei modulations (0) is plotted vs. the mass of the perturber. To 
show the effect on the shape of the resulting ei modulation, three extreme cases 
are presented in the bottom panels: 1) M3 = 2Mj and P3/P2 ~ 23; 2) the solution 
with maximum eccentricity amplitude, with M3 = 48Mj and P3/P2 = 13.4; 3) 
M 3 = 152Mj and P3/P2 = 36. right) The light curve of FS Aur folded with the 875 
day period (top panel). In the bottom panel, the eccentricity modulation is shown 
for the best solution to explain the VLPP (M 3 = 48Mj and P3/P2 = 13.4). 



of FS Aur on comparable time scales. The extended description of numerical 
calculations and their results and implications upon observed peculiarities of FS 
Aur are presented in Chavez et al. (2010) [7] . 

3. A precessing, magnetic WD in FS Aur 

Tovmassian ct al. (2007) proposed that FS Aur possesses a rapidly rotating 
magnetic WD precessing with the LPPp]. Neustroev et al. (2005) intended to 
reveal the spin period of the WD by fast optical photometry [S]. They found a ~ 
100 sec peak in the power spectrum, but it was not conclusive evidence after taking 
into account the weakness of the signal and the peculiarity that it was clearly 
visible only half of the OP and practically absent in the other half. Nevertheless, 
Tovmassian et al (2007) discovered the LSP[?] and found no better way to explain 
it than the same model of the precessing, magnetic WD, which creates a bright 
source of emission at the inner edge of the truncated accretion disc around the spot 
where the coupling of the disc matter with the magnctosphcre of the WD occurs. 
We conducted X-ray observations using Chandra and Swift X-ray telescopes to 
detect the magnetic WD, determine its spin period and check the proposed model. 
Here, we report preliminary results of these observations. FS Aur is a relatively 
bright X-ray source with a count-rate of 0.38 cnts/s as observed with Chandra's 
ACIS-S3 (back-illuminated) CCD array during a 25 ks exposure. The spectrum is 
similar to a CV and is best fitted with the multi-temperature mkcflow model with 
a range of temperatures from 0.3 to 3.8 keV. The spectrum is highly absorbed 
with the best fit giving a hydrogen column density of Nh = 0.4 x 10 22 cm~ 2 . 
The non-magnetic CVs usually have absorption column density of a few times 
10 20 cm -2 [9]. A highly absorbed spectrum is a good indication of a Intermediate 
Polar (IP) nature for the object[10]. However, the X-ray LC does not show any 
clear spin period spike in the power spectrum. There is nevertheless a marginal 
peak at 101 sec, which is the presumed spin period of WD in FS Aur, and it 
corresponds to a similarly weak signal in the optical LCs, which might be an 
indication of the reality of the peak. The weakness or absence of a spin period 
signal in X-rays may sound worrisome for claims that FS Aur is an IP and contains 
a moderately magnetic WD, but Ramsay et al (2008) showed that in the case where 
the magnetic pole and spin axis are aligned within a few degrees, it will give rise 
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Fig. 4. The X-ray light curve in soft and hard energy channels. 



to an X-ray modulation less th an a few percent, which would be undetectable 
due to intrinsic X-ray flickering [ 1 Oj . Actually, the model involving a precessing, 
magnetic WD in order to explain the LSP would work if the magnetic pole of the 
WD is close to the spin axis of the star. 

The really important result of the X-ray observations is the detection of OP 
modulation of the X-ray LC. As a matter of fact, the X-ray modulation has a 
pulse profile similar to IPs (i.e., FO Aqr QT|). OP modulation is detected in al- 
most 50% of IPs [12] with an amplitude of modulation decreasing towards higher 
energies, and the modulation of FS Aur is stronger in the soft band, indicating 
photoelectric absorption typical of IPs. In Fig 07] the X-ray LC of FS Aur is pre- 
sented in two energy bands: soft (0.3-2.0 keV) and and hard (2-10 keV). What 
is remarkable is that one out of five pulses seems to be missing in the LC. The 
Baycsian analysis of photon distribution shows that the absence is statistically 
significant. In fact, XY Ari exhibits X-ray modulation at the OP that tends to 
appear and disappear, which they explain with precession, but of tilted/warped 
disc[T3]. We believe in case of FS Aur, it is rather evidence supporting the hy- 
pothesis of WD precession. Of course, one missing pulse is not good evidence, 
but follow up X-ray observations with the Swift telescope provided more signifi- 
cant evidence in favor of WD precession. The Swift observations provided shorter 
continuous segments of the LC, but a larger amount of total coverage. They are 
comprised of short runs of about 700 — 1 200sec each every ~ 6 000 sec with total 
coverage of about 200 ks in three large blocks. The sensitivity of the telescope 
is lower than that of Chandra (the object is detected as a 0.09 cuts/ sec source), 
but the resulting spectrum and LC are remarkably similar, except the LC has 
a marked phase shift compared to the LC obtained by Chandra! In FiglBT] the 
X-ray LC of FS Aur folded with the OP is presented. The upper panel presents 
Chandra data and the lower is Swift data, they are both binned into equal time 
bins with phases calculated from the same zero point determined from the optical 
spectroscopy. The shift between the two curves is about 0.2 orbital phases. Nei- 
ther zero point nor the amount of shift have any special relevance; however, the 
presence of the shift is so far the strongest evidence in favor of precession of the 
WD in FS Aur. 

4. Conclusions 

FS Aur is a 85.797 min DN, which, in addition to the orbital period, demon- 
strates a number of longer than orbital spectroscopic and photometric periods. 
We demonstrate here that in addition to LPP and LSP, the object also shows a 
very long photometric period, varying up to 2 magnitudes around the quiescent 
average magnitude with a period of ~ 900 days. Such a long period may be ex- 
plained by the presence of sub-stellar third body on a circular orbit around the 
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Fig. 5. left) The power spectrum of X-ray light curve, right) The light curves from 
Chandra and Swift folded with the OP and the same To. 



CB with orbital period much smaller than the detected VLPP. The presence of a 
third body might be helpful in understanding not only the VLPP but also other 
unusual features in the LC of FS Aur on shorter time scales since the eccentricity 
modulation occurs on three different time scales. In the case that the mass of 
third body is somewhere in between 25 and 65 Mj, the amplitude of the VLPP 
will be much larger than on short time scales. Among other things, it is probably 
be more natural to expect a precessing white dwarf in this particular CV, although 
it is difficult to quantify influence of the third body on a rapidly rotating WD. 

The precession hypothesis, proposed by Tovmassian et al. (2007) finds addi- 
tional support from X-ray observations of FS Aur. The hard X-ray spectrum and 
the form of the modulated X-ray LC by the orbital period provide evidence for 
the suggestion that FS Aur is an IP. In IPs, the orbital modulation is believed 
to be a result of phase-dependent photoelectric absorption. If we consider that 
the coupling region of the accretion disk material with the magnetosphere of the 
WD, where the absorption of X-rays takes place, is the same spot in the inner 
edge of accretion disk that produces the high RV wings of the emission lines in 
optical, then the precession hypothesis makes a lot of sense. Since that spot will 
go around following the precessing magnetic pole of the WD, it will cause phase 
shifts in the X-ray LC. Thus, we think that the detection of these shifts provides 
excellent evidence for precession and indirectly confirms that the weak ~ 100 sec 
peak in power spectra detected in optical and X-rays is indeed the spin period of 
the WD. FS Aur has proved to be very complex, unusual even for a CV, case. It 
seems that a combination of a triple system containing an IP might resolve most 
of the issues, even though some mechanisms are not readily understood. 
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